Dewetting of thin films on heterogeneous substrates: Pinning vs. coarsening 
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We study a model for a thin liquid film dewetting from a periodic heterogeneous substrate (template). The 
amplitude and periodicity of a striped template heterogeneity necessary to obtain a stable periodic stripe pattern, 
i. e. pinning, are computed. This requires a stabilization of the longitudinal and transversal modes driving the 
typical coarsening dynamics during dewetting of a thin film on a homogeneous substrate. If the heterogeneity 
has a larger spatial period than the critical dewetting mode, weak heterogeneities are sufficient for pinning. A 
large region of coexistence between coarsening dynamics and pinning is found. 

PACS numbers: 68.15.4e, 68.55.-a, 47.20.Ma, 47.20.Ky 



Templating and controlled rupture of liquid films on chem- 
ically structured substrates have provoked many experimental 
efforts [Q, ^ pi 0, BL ^] but the conditions for the desired imag- 
ing of the template structure onto the deposit pose many open 
questions. First theoretical results for liquid layers on struc- 
tured surfaces with strong heterogeneities of wettability, i.e., 
stepwise alternating hydrophilic and hydrophobic stripes in an 
aqueous system, describe different morphological transitions 
when changing the size of the heterogeneous patches or the 
volume of the deposited liquid [^, |^, ^ [l0| . Experimentally, 
the challenge of preventing the dewetting pattern from coars- 
ening has been met by evaporation of the solvent JTl|], by lit- 
erally freezing the system [|l2[ or by using a heterogeneous 
substrate 

In the present paper we study the transition between coars- 
ening and pinning for thin films on weakly structured sub- 
strates that possess spatial modulations of the molecular in- 
teraction terms. Thin films on homogeneous substrates can 
be unstable to spinodal dewetting, see e.g. [|l3[ [l^j. Sharma 
and coworkers proposed a model that contains polar and apo- 
lar components of molecular interactions and reproduced the 
dynamics of dewetting thin films jl5|]. Recently, it has been 
shown that this model has periodic stripe solutions that are 
unstable to coarsening flip]]. Here we use a slightly differ- 
ent model derived from diffuse interface theory by Pismen 
and Pomeau jl7| ] as a starting point. We emphasize that this 
model has very similar dynamics [ |l8| ] and stationary solu- 
tions fll^ ] as the Sharma model. The impact of a heteroge- 
neous substrate on the stationary film profiles, their linear sta- 
bility and resulting conditions for morphological transitions 
are presented. The dynamics of thin films on periodic stripe 
templates have been studied in JlOt l with the original Sharma 
model by numerical integration and for a given heterogeneity 
of large amplitude. Using bifurcation and stability analysis, 
we can provide, however, a more systematic study of the ef- 
fects of parameters like average film thickness h, amplitude e 
and spatial period Phet of the heterogeneity. 

Our starting point is a homogeneous or weakly heteroge- 
neous substrate. We choose parameters in the unstable regime 
(cross in Fig. [jja)). Then a striped solution of a spatial pe- 
riod P = P^t is considered. For illustration we display a 
system of length 2Ph e t (see Fig. |l|(b)-(e)). The periodic stripe 
pattern (Fig. [0(b)) is now unstable against several transversal 
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FIG. 1 : (a ) Phase diagram for thin films on a homogeneous substrate 
(after [|19|]). The cross shows the parameter values for which we 
present results in detail, (b) Schematic display of the template (dot- 
dashed line), the periodic film profile with the same spatial period 
(solid line). Shaded areas show where h > h. (c) Initial stage of 
the transversal instability and (d),(e) final stages of variants of the 
longitudinal coarsening instability. 

(Fig. [l|(c)) and longitudinal (Figs. [j](d),(e)) perturbations. Sta- 
bility analysis allows us to trace the corresponding eigenval- 
ues. Increasing the strength of the heterogeneity, they all will 
acquire negative real parts and the desired pattern in Fig. [jjb) 
is pinned, i.e. becomes stable against small perturbations. 

The evolution equation for the film thickness profile 
h(x, y, t) contains a disjoining pressure that is derived by 
combining the Stokes equation in long wave approximation 
Q2Q] with diffuse interface theory [|l7||: 
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with a free energy f(h, x) — k(x) e (e — 2) 
and the ratio G of gravitation to mean molecular interactions 
as well as the spatially varying strength of the molecular in- 
teractions k(x), a > is a small parameter describing the 
wetting properties in the regime of partial wetting [ |l7[ [lf-fl . 
This dimensionless form is obtained after scaling the origi- 
nal quantities as in Ref. [|l9|]. The molecular interactions that 
become important on the nanometer scale [^TJ] are incorpo- 
rated through the disjoining pressure H(h) contained in the 
derivative of the free energy dhf(h, x) — —n{x)Ii{h) + Gh. 
The chosen H(h) is qualitatively equivalent to other pressures 
consisting of destabilizing short-range and stabilizing long- 
range interactions [ p0[ | but does not suffer from the usual di- 
vergency for vanishing film thickness JT^]. In the absence 
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of heterogeneity, the model possesses two control parameters, 
the dimensionless quantity G and the average film thickness 
h that is chosen as conserved quantity. A phase diagram in G 
and h indicating the region of spinodal dewetting is shown in 
Fig. 0(a). 

Heterogeneous substrates with a smooth change in the wet- 
tability are modelled by a spatial sinusoidal modulation of the 
overall strength of the disjoining pressure 



k(x) = 1 + ecos(2irx /Phet) 



(2) 



The system behaviour is controlled by the amplitude, e, and 
the imposed periodicity, Phet of the heterogeneity. Locally 
small k correspond to effectively larger G [Jl9|]. Droplets are 
preferably located at minima of k(x). 

We choose the system size, L, as a multiple of Phet, apply 
periodic boundary conditions and vary e > 0. The Lyapunov 
functional 



E = - 



f(h,x) 



dx 



(3) 



gives the energy of a striped solution h(x,t) [|l9|. During 
the time evolution of a given initial film thickness profile this 
energy decreases and eventually settles in a minimum when 
the system approaches a linearly stable stationary solution of 
Eq. ([y). One can_determine these solutions directly by setting 



dth = in Eq. ([!]) and integrating twice [Q, yielding 

= d xx h - d h f(h, x) + Ci 



(4) 



The integration constant C\ and the spatial period P now 
parametrize a two-parameter family of periodic solutions, 
i.e. periodic patterns of holes and droplets, that can be cal- 
culated by continuation techniques p3|]. P can be equal 
to Phet or its multiples. Focusing on situations with con- 
served liquid volume we adopt the mean film thickness h = 
1/L Jq h(x)dx as the second parameter beside the period, 
using C\ as a Lagrange multiplier. Typical film profiles are 
shown in Fig. ||(a),(b). 

For homogeneous substrates with e = 0, flat films are 
stationary solutions that, depending on mean film thickness 
and the parameter G, may be unstable to infinitely small 
(spinodal dewetting) or finite (nucleation) disturbances of the 
film surface. Rupture due to the spinodal mechanism oc- 
curs for perturbations with periods larger than a critical value, 
P c = 2ir/k c with k\ = —dhhf(h,x)\h [p^|]. On a slower time 
scale the initially formed holes coalesce, the pattern becomes 
coarser and tends to the absolute minimum of the energy at 
the largest possible P, i.e. the system size. 

( i) Perturbation theory - If the heterogeneity is switched on, 
the flat film is no longer a solution to Eq. (Q). However, for 
very small e< lan analytical expression for the weakly vary- 
ing stationary film profile can be calculated. Rewriting the 
heterogeneity as e(e ikh " tX + c.c.) jl with khet = 2ir/Phet and 
using the volume conserving ansatz h(x) = h + a(e lkh " tX + 
c.c.)/2 with 0(a) — 0(e) in Eq. (||), gives to first order in e 
a stationary solution with 
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FIG. 2: (a) Stationary film profiles to Eq. [j] with P = P het (thick) 
and P = 2Ph e t (thin), (b) Portraits of all solutions with P — Phet- 
Only the solid profile is linearly stable, (c) The heterogeneity and (d) 
the two coarsening modes due to translation (dashed) and transfer 
of mass (dotted) as sketched in (a). Parameters are e = 0.1, G = 
0.1, h = 2.5, P h et =50. 

This solution is only valid for a<l, For P c < Phet, the so- 
lution is modulated in phase with the heterogeneity, whereas 
for P c > Phet the phase is shifted by ir. Results obtained by 
perturbation theory are depicted in Fig. |](b). 

( ii) Numerical bifurcation results - Suppose a variety of het- 
erogeneous substrates with appropriate period Phet but dif- 
ferent e is available and a dewetting pattern with spatial pe- 
riod P = Phet is desired. We choose a linearly unstable 
mean film thickness h such that the critical wavelength of 
spinodal dewetting P c is of the same order as Phet- As an 
example, we choose P c « P/j et /1.5 with Phet — 50 and 
G = 0.1, h = 2.5. This ratio of Phet/Pc is known to give 
good templating for large heterogeneities JTo| ] . The stationary 
solutions with P = Phet and P = 2Phet are computed [ p2[ ] 
as e is increased using the weakly modulated solutions as a 
starting point for the continuation [H]. Figs. ||(a),(b) and||(a) 
showpro files of the solutions and the bifurcation diagram and 
Fig. 0(b) compares the results with the small e approximation 
Eq. (Eh. The weakly modulated solutions are indeed in line 
with the results of perturbation theory. Breaking the transla- 
tional symmetry at e = gives rise to two branches of solu- 
tions starting at the large amplitude solution known from the 
homogeneous case [p9[]. Only the branch of largest amplitude 
possesses a phase shift. The branches of solutions in phase 
cease to exist in a saddle-node bifurcation at e « 0.22. 

(Hi) Energy and longitudinal stability - Thick curves in 
Fig.|(a) compare the energies (calculated with Eq. (||)) of the 
different solution branches given in Fig. |3|(b). For the chosen 
parameter values the solutions in phase have always a larger 
energy than the ones out of phase. When increasing e the 
energy of the lower branch decreases further indicating that 
the heterogeneity favors this pattern with P = Phet- Thin 
curves in Fig. |3|(a) denote solutions with P = 2P/ let . To gain 
more insight into the conditions of pinning the pattern by a 
heterogeneity, we analyse the linear stability of the stationary 
solutions in systems of different sizes [f24|]. For the stability 
analysis of periodic solutions, one usually employs a Floquet 
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FIG. 3: (a) Relative energy of stationary solutions to Eq. (|l|) with 
P = Phet (thick) and P = 2P het (thin) versus e. A flat film 
at e — yields » 0.155. Solid curves correspond to coex- 

isting linearly stable solutions with, (b) Bifurcation diagram rep- 
resenting maximum (upper half) and minimum (lower half) of the 
solutions with P = Phet- The approximation (^|) is denoted by 
the thin dot-dashed lines and other line styles match with (a) and 
Fig. ^(b). (c) Eigenvalues A with largest real part of the lowest en- 
ergy solutions with P — Phet (thick solid branch in (a)). Solid 
curves correspond to Goldstone modes at e = and broken curves 
to interaction modes with period 2P. Line styles match the pro- 
files in Fig. ^d). Triangles correspond to period doubling bifur- 
cations where solutions with P — 2Phet emerge. Parameters are 
a = 0.1, G = 0.1, h = 2.5, P het = 50 which yields a/e » 0.7. 

ansatz. This enables us to get the stability of periodic solutions 
in large systems corresponding to large ratios n = LjPhet- 
We find however that the most dangerous, potentially unstable 
longitudinal modes that induce coarsening of dewetting films 
are already present in systems for n = 2. Hence, we constrain 
the discussion to system sizes L < 2Pf let in x— direction. In 
the transversal y-direction, we had to use much larger length 
w 30Phet for a stability analysis, see below. 

In a system of size equal to one period of the heterogeneity, 
i.e. L = Phet, the entire lower branch in Fig. ^|(a) is lin- 
early stable, whereas the two other branches are unstable. For 
e = the linearly stable solution possesses two zero eigenval- 
ues (Goldstone modes) corresponding to translation symmetry 
and mass conservation. For c ^ the translation symmetry 
is broken and the corresponding eigenvalue becomes negative 
(thin solid line in Fig. j^c)). Mass conservation is maintained 
and the corresponding eigenvalue remains zero for all e (thick 
solid line in Fig. ||(c)). 

In a larger system with L = 2Pf let , the first stage of coars- 
ening can be studied which is the fastest since nearest neigh- 
bours interact and the energy gain is largest. If templating 
shall be successful then this first stage needs to be prevented 
by pinning. In the stability analysis of the two periods, two 
new eigenvalues appear that correspond to asymmetric com- 
binations of the Goldstone modes (see Fig. 0(d)). They rep- 
resent two possible modes of coarsening: (a) Shift of droplets 



towards each other caused by a combination of opposite trans- 
lational modes (dashed arrows in Fig. ||(a) and dashed curves 
in Figs. ||(d) and||(c)); (b) Mass transfer between neighbour- 
ing unmoved droplets caused by a combination of opposite 
volume modes (dotted arrows in Fig. ||(a) and dotted curves in 
Figs. H^d) and |](c)). These modes have first been recognized 
in the Cahn-Hilliard equation [p6|]. For small e both eigenval- 
ues are positive, implying the instability of the wanted pattern 
to coarsening. The mass transfer proceeds faster than the shift 
of droplets. But, as e increases, both become negative, ren- 
dering the solution with period Phet linearly stable. At the 
two crossing points period doubling bifurcations occur where 
stationary solutions ofperiod P = 2Ph e t emerge (compare 
thin branches in Fig. ||(a)). Altogether, the desired pinning 
solution with P = Phet is linearly stable for e > 0.002. 

The value of this critical e depends on parameters and the 
specific form of H(h). Since the two symmetries connected 
with the coarsening modes are present for arbitrary choices of 
H(h) the result does not depend on this choice. Accordingly, 
two pairs of solutions with P = 2Ph e t appear in any given 
system at finite e. 

These four branches of solutions with P = 2Ph e t (thin 
curves in Fig. |3|a)) have the following linear stability. Solu- 
tions on the dotted branch carry two positive eigenvalues. The 
short-dashed branch has still one positive eigenvalue that leads 
to a shift of the pattern towards solutions of the solid branch. 
The long-dashed branch has one positive eigenvalue and is a 
saddle that divides evolutions by translation of two droplets 
towards one on the solid branch or back to the P = Phet 
branch. The entire thin solid branch is linearly stable in 
L = 2Phet and represents the coarse solution competing with 
the desired pattern (see Fig. [j](e) resp. the thin profile in 
Fig. ||(a)). Here, at large h, the coarse profile of lowest en- 
ergy emerges from the translation mode, whereas for smaller 
mean film thicknesses (h < 2 for G = 0.1) the coarse profile 
corresponding to the transfer mode is stable (see Fig. [l|(d)) 
and has lowest energy [p5||. 

(iv) Coarsening versus pinning - As seen in (Hi) coarsening 
is favored for e < 0.002 and solutions with P = 2P/ let are 
the only stable ones. Consequently they also have the lowest 
energy (for L = 2P/ let ). For e > 0.002 multistability be- 
tween linearly stable solutions with P = Phet and P = 2Ph e t 
occurs and the initial condition becomes crucial for the final 
dewetting structure. Below e w 0.157 the coarse solution has 
lowest energy, whereas at larger e the desired pinned pattern 
is energetically favored. Choosing suitable initial conditions 
enables pinning at much smaller e. But coarsening will still 
occur at larger e if initial conditions are chosen accordingly. 
At even larger e > 0.5 the pinned pattern is the only possibil- 
ity because the coarse solutions cease to exist. 

Summarizing the results for a range of values Phet at fixed 
G = 0.1 and h = 2.5 we obtain the "morphological phase di- 
agram" Fig. |j. The product e * Phet emerges as a rough funda- 
mental parameter: Coarsening prevails when it has low values 
while for large values the pattern pins to the heterogeneity as 
desired. At intermediate values multistability is found where 
the initial condition selects the final outcome. 

(v) Transversal stability - The transversal stability of the 
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FIG. 4: Morphological phase diagram of templating showing regions 
in the parameter plane (e, Phet) with different behavior of the thin 
film on a heterogeneous substrate. The shaded band separates pa- 
rameters of ultimate coarsening from ultimate pinning. Inside the 
shaded band multistability is present with the desired pattern being 
the energetically minimum inside the dark shaded area. Parameters 
are a = 0.1, G = 0.1, h = 2.5, P c ~ 33. The triangles correspond 
to the equivalent symbols in Fig. |3](a). 

pinned profiles (P = Phet) can be obtained from the eigenval- 
ues A of small perturbations Sh(x) exp(Xt)(exp(iky) + c.c.) . 
We observe a long wavelength instability at small e on a length 
scale much larger than the longitudinal modes. The most dan- 
gerous branch of transversal perturbations emerges from the 
zero-eigenvalue representing mass conservation in the film. 
For e = 0.001 for instance, we found that the fastest growing 



mode had a wavelength of ks 900 = 18Phet (see Fig. [j](c)). 
The transversal modes are stabilized for parameters above the 
solid curve in Fig. || For the present choice of 11(h) and G, h 
the transversal instability restricts the stability range of the 
pinned solution more than longitudinal coarsening. 

(vi) Summary - We have studied the conditions for success- 
ful templating or pinning of thin films that are unstable to spin- 
odal dewetting. The main result is that a periodic stripe pattern 
can be pinned if the heterogeneous substrate suppresses the 
transverse and longitudinal instabilities experienced by such 
patterns on a homogeneous substrate. While previous stud- 
ies have employed direct numerical simulation to study the 
dynamics of thin films, we have used numerical bifurcation 
and stability computations. This allows an efficient scanning 
of parameters characterizing the heterogeneous substrate and 
the film dynamics. Comparing the two length scales P c and 
Phet of spinodal dewetting and heterogeneity, respectively, we 
find that patterns are not pinned to heterogeneities on a much 
smaller scale than P c , and that smaller P c need weaker het- 
erogeneities to pin the pattern. In consequence, templating 
can be best controlled by choosing the initial mass of fluid 
which yields the smallest P c , i.e. the film thickness where 
the derivative —dhhf(h, x)|^ is maximal (for constant G), see 
also [[27]]. The transition from coarsening to pinning is hys- 
teretic, giving rise to an extended region where the desired 
pinned pattern coexists with coarser structures. Our results 
and the methodology also apply to other strategies of stabiliz- 
ing a periodic pattern including introduction of anisotropy or 
convective flows [E8|] ■ 



* Electronic address: 
' Electronic address: 



baer@mpipks-dresden.mpg.de 



thiele@mpipks-dresden.mpg.de 

M. Yoon, and 



[1] L. Rockford, Y. Liu, P. Mansky, T. P. Russell, 

S. G. J. Mochrie, Phys. Rev. Lett. 82, 2602 (1999). 
[2] M. Mertig, R. Kirsch, W. Pompe, and H. Engelhardt, Eur. Phys. 

J. D 9, 45 (1999). 
[3] H. Gau, S. Herminghaus, P. Lenz, and R. Lipowsky, Science 

283, 46 (1999). 
[4] D. E. Kataoka and S. M. Troian, Nature 402, 173 (1999). 
[5] M. Gleiche, L. F. Chi, and H. Fuchs, Nature 403, 794 (2000). 
[6] A. M. Higgins and R. A. L. Jones, Nature 404, 476 (2000). 
[7] P. Lenz and R. Lipowsky, Phys. Rev. Lett. 80, 1920 (1998). 
[8] C. Bauer and S. Dietrich, Phys. Rev. E 61, 1664 (2000). 
[9] R. Konnur, K. Kargupta, and A. Sharma, Phys. Rev. Lett. 84, 

931 (2000). 

[10] K. Kargupta and A. Sharma, Phys. Rev. Lett. 86, 4536 (2001). 
[11] U. Thiele, M. Mertig, and W. Pompe, Phys. Rev. Lett. 80, 2869 
(1998). 

[12] J. Bischof, D. Scherer, S. Herminghaus, and P. Leiderer, Phys. 

Rev. Lett. 77, 1536(1996). 
[13] G. Reiter, Phys. Rev. Lett. 68, 75 (1992). 
[14] S. Herminghaus etal, Science 282, 916 (1998). 
[15] A. Sharma and R. Khanna, Phys. Rev. Lett. 81, 3463 (1998). 
[16] U. Thiele, M. G. Velarde, and K. Neuffer, Phys. Rev. Lett. 87, 

016104 (2001). 

[17] L. M. Pismen and Y. Pomeau, Phys. Rev. E 62, 2480 (2000). 



[18] M. Bestehorn and K. Neuffer, Phys. Rev. Lett. 87, 046101 
(2001). 

[19] U. Thiele, M. G. Velarde, K. Neuffer, and Y. Pomeau, Phys. 

Rev. E 64, 031602 (2001). 
[20] A. Oron, S. H. Davis, and S. G. Bankoff, Rev. Mod. Phys. 69, 

931 (1997). 

[21] J. N. Israelachvili, lntermolecular and Surface Forces (Aca- 
demic Press, London, 1992). 

[22] We computed all stationary solutions with P — nPhet and 
n < 10 but the qualitative picture is not changed. 

[23] E. Doedel, A. Champneys, T. Fairfrieve, Y. Kuznetsov, B. Sand- 
stede, and X. Wang, AUT097: Continuation and bifurcation 
software for ordinary differential equations (Concordia Univer- 
sity, Montreal, 1997). 

[24] Space has been discretized by 64 Fourier modes and the spec- 
trum of eigenvalues as well as the corresponding eigenmodes 
are calculated from the linearization of Eq. (jl]) around the sta- 
tionary solutions. 

[25] Lutz Brusch, Heiko Kiihne, Uwe Thiele, and Markus Bar, in 
preparation. 

[26] K. Kawasaki and T. Ohta, Physica A 116, 573 (1982). 
[27] K. Kargupta, R. Konnur, and A. Sharma, Langmuir 16(26), 
10243 (2000). 

[28] A. A. Golovin, A. A. Nepomnyashchy, S. H. Davis, and 
M. A. Zaks, Phys. Rev. Lett. 86, 1550 (2001). 



